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Abstract 

We have studied holographic superconductors with spherical symmetry in the 
Horava-Lifshitz gravity by using a semi analytical method, and also we have cal- 
culated the critical temperature and shown when the condensation will appear in 
a similar pattern as in the Einstein-Gauss- Bonnet gravity. We have computed the 
dependency of the conductivity as a function of frequency in this new non-relativistic 
model of quantum gravity. 
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I. INTRODUCTION 

As a phenomenological fact, superconductivity is usually modeled by a Landau- Ginzburg 
Lagrangian where a complex scalar field develops a condensation in a superconductive phase. 
To have a scalar condensation in the boundary theory, Horowitz and his collaborators [l| 
introduced a U(l) gauge field and a conformally coupled charged complex scalar field in 
the black hole background. That potential corresponding to the conformal mass is nega- 
tive, although above the Breitenlohner-Freedman (BF) bound [2i] it does not cause any in- 
stability in the theory. To solve the negative mass problem Basu et.al and his collaborators 
jsj showed that the presence of the vector potential effectively modifies the mass term of the 
scalar field as we move along the radial direction r and allows the possibility of developing 
hairs for the black hole in some parts of the parameter space. In their model there was no 
explicit specification of the Landau- Ginzburg potential for the complex scalar field. The 
development of condensations relies on a more subtle mechanism violating the no hair theo- 
rem. Further Wen investigated the holographically dual description of superconductors in (2 
+ l)-space time dimensions in the presence of inhomogeneous magnetic field and observed 
that there exist type I and tjnpe II superconductor [4]. the existence of holographic super 
conductors was established in[l|, |5[. From the (d dimensional) field theory point of view, su- 
per conductivity is characterized by condensation of a generally composite charged operator 
O in low temperatures T < .In the gravitationally dual (d+1 dimensional ) description 
of the system, the transition to the super conductivity is observed as a classical instability 
of a black hole in an anti-de Sitter (AdS) space against perturbations by a charged scalar 
field ip. The instability appears when the black hole has Hawking temperature T = Tc. For 
lower temperatures the gravitational dual is a black hole with a non vanishing profile for the 
scalar field ip. The AdS/CFT correspondence relates the quantum dynamics of the boundary 



operator O to a sinaple classical dynamics of the bulk scalar field ip js, 7|. Following Hart- 
noil et al works in and also Maeda and Okamura [Sj, we will find out they studied the 
perturbation of the gravitational system near the critical temperature Tc, and they obtained 
the superconductor's coherence length via AdS/CFT (antide Sitter/conformal field theory) 
correspondence, and also they added a small external homogeneous magnetic field to the 
system, and found a stationary diamagnetic current proportional to the square of the order 
parameter being induced by the magnetic field. Their results agree with Ginzburg-Landau 
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theory and strongly support the idea that a superconductor can be described by a charged 
scalar field on a black hole via AdS/CFT duality. From a pure classical treatment, there 
is more efforts to deal with BH in Ads backgrounds. Black holes in anti-de Sitter (AdS) 
spacetime in several dimensions have been recently studied. One of the reasons for this 
intense study is the AdS/CFT conjecture stating that there is a correspondence between 
string theory in AdS spacetime and a conformal field theory (CFT) on the boundary of that 
space. For instance, the M-theory on AdS"^ x S''' is dual to a non-Abelian super conformal 
field theory in three dimensions, and type IIB superstring theory on AdS^ x seems to be 
equivalent to a super YangMills theory in four dimensions ^. 
Recently, a power- counting 



renormalizable, ultra-violet (UV) complete theory of gravity was 
proposed by Horava in lCHl3l|. Although presenting an infrared (IR) fixed point, namely 
General Relativity, in the UV the theory possesses a fixed point with an anisotropic, Lif- 
shitz scaling between time and space of the form — )■ £ x*, t — )■ t, where i, z, and t 
are the scaling factor, dynamical critical exponent, spatial coordinates and temporal coor- 
dinate, respectively. According to the Bias et al arguments [ij], it seems that this model 
must be modified by some terms to avoid from strong coupling, instabilities, dynamical in- 
consistencies and unphysical extra mode. As we know that there are two explicit families 
of exact solutions for a spherically symmetric background without projectability condition 
in HL gravity and other solutions all are the familiar GR solutions i.e Ads^-Schwarzschild 
solutions. First solution belongs to the |l5| known asymptotically fiat KS solution and as 



we have shown that in spite of the GR BHs, its timelike geodesies is stab 



The other 



non trivial solution was found by Lu-Mei et al. [17|], and recently Tang 18|] investigated the 
general solutions of the HL theory under both projectability and non projectability condi- 
tions. His paper contains all the former solutions and at the end of it, he presented two 
new families of exact solutions - only in a neutral case- which both of them are valid in the 
corner of the validity of the IR limit of the HL theory i.e A = 1 and these solutions can be 
interpreted as two new forms of the BHs for HL gravity. 

Recently the works were done about the Holographic Superconductors for a new topo- 
logical BH in HL gravity describing a topological black hole solution whose horizon has an 



arbitrary constant scalar curvature 



19 



20 



22[ |. They found that it is more applicable for 



the scalar hair forming, when the parameter of the detailed balance( e) becomes larger, and 
harder when the mass of the scalar field is larger. Also they calculated the ratio of the gap 
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frequency in the conductivity with respect to the critical temperature. Briefly they inves- 
tigated the effects of the mass of the scalar field and the parameter of the detailed balance 
on the scalar condensation, the electrical conductivity, and the ratio of the gap frequency in 
the conductivity at the critical temperature. 

There are many interesting features for critical phenomena and superconductivity when we 



are working on hig^ 
Bonnet corrections 
of works 



ler orders corrections, specially when we are interesting in the Gauss- 



2l[ |. The same phenomenology has been discussed by Wang in series 



20l |. These phenomena and it's physical consequences are very similar with our 



analysis in the HL theory and we can generalize their results to our higher order theory in 
the non relativistic regime. 

In this work we have discussed a type of solutions which has been reported inp/7j. In 
Sec. 2 we have presented spherically symmetric black holes' solutions in Horava- Lifshitz 
gravity with the action without the condition of the detailed balance. In Sec. 3 we have 
explored the scalar condensation in the Hof ava-Lifshitz black hole by analytical approaches. 
In Sec. 4 the matching solutions and the critical temperature have been found. In Sec. 5 
we have computed the conductivity of our model and shown the behavior of the real part 
of the conductivity as a function of frequency per tempereture. We have summarized and 
discussed our conclusions in the last section. 



II. SOLUTIONS OF THE HOrAVA- LIFSHITZ GRAVITY 

Since in the HL theory, the dynamical quantities are the shift Ni{t,x), lapse N{t,x) and 



metric hij] therefore in the ADM formalism 23| : 



ds^ = -N^df + hij{dx' + N'dt){dx^ + NHt) (1) 

If we restricted ourselves to the static metrics hij, there is two possibility for the time 
dependency of the two remaining functions. In many cases as in Lu-Mei-Pope 171], we can 
relax the shift function by a formal going to the Schwarzschild gauge and rewriting the 
static solution with spherical symmetry in GR. Thus for solutions in the usual Schwarzschild 
gauge the only function is the lapse. According to the terminology of the Horava theory, a 
projectable solution is a solution with a time dependent lapse and a non projectable one is a 
vise versa. Many authors consider the non projectable version as an exact solution. Another 
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problem returns to the choice of the potential term. The first choice is due to the detailed 



balance principle 



24j |. but in the original work of the Horava in the context of the cosmology 



this principle implies a negative cosmological constant in contrary with the observational 
evidences. The other problem is avoiding from the ghost excitations [14] restricting one 
to accept a value of the A < | or A > 1. Instability and strong coupling impose another 
difficulties for it. Far from all of these problems we rewrite an explicit spherical symmetric 



solution for HL theory following Lu-Mei-Pope work 



17|. 



A. New static neutral BH solution 



Following the ADM formalism, the action of this HL gravity with a soft violation of the 
detailed balance condition is given by: 



S = dtd^x^N{CK - 
Jm 



(2) 



K 



4(3A - 1) 



W) 



8(3A - l)R 



The at are the coupling parameters [17| , and Cij is the Cotton tensor [12| . With the metric 
ansatz as in 



ds^ = -Nirfdt^ + -rr^idr + N^'dtf + r^dU'^ 



The following solution in the UV region has been found 17|: 



AT'- = 
2A ± V6A - 2 
A^l 
7 = 5-1 



(3) 



(4) 



fir) = / = 1 - 



A 



W 2 5 

r — ar 



N{r) = N = fir"'^^f 



(5) 
(6) 



where are constants. This solution is asymptotically AdS^" and thus it is useful in 
the AdS/CFT correspondence scenario for the Holographic superconductivity. The Hawk- 



ing temperature is given by the usual Gibbons-Hawking calculus 



25| . therefore the Unruh 
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temperature can be written in the form 26| : 



T 



27r 



\r=rH 



An 



h-y'{h) 



A 

47r 



(7) 



in order to satisfy the positivity of the temperature, we must require f3 < when both Aw 
and a are positive simultaneously. 

III. FIELD EQUATIONS FOR SCALAR CONDENSATION SCENARIO 

FoUowingthe HartnoU, Herzog and Horwitz general framework to the holographic super- 
conductors [1, in the limit where the scalar field does not back- react on the geometry 



the solution for the background geometry is that of the dyonic black hole j28| . In this paper, 



the charge density of the background 



13, 



271] is neutral, so both the electric and mag- 



netic charge of the dyonic black hole have been set to zero. The Maxwell-scalar sector is 
decoupled from the gravity sector, therefore the minimal ingredients we need to describe a 
holographic superconductor are conserved energy momentum T'^", Global U(l) symmetry, 
conserved current J'^ and finally charged operator O condensing at low temperature (/i, u 
runs over t, x, y). The most basic entries in the AdS/CFT dictionary js, 7| tell us that there 
is a mapping between field theory operators and fields in the bulk . In particular, T^^ will 
be dual to the bulk metric Qab-, the current will be dual to a Maxwell field in the bulk 
Aa, and the dual of charged scalar field ip is O { here a, h runs over t, x, r). We can now 
study the Maxwell-scalar theory in the black hole background with Lagrangian: 



-\F^-\^^-^A^\^ + 2i^ 



(8) 



The only dimensional parameter in the Lagrangian is L related to the AdS radius, and the 
full set of equations of motion for the fields ip and are : 



1 



'-9 








(9) 
(10) 



respectively, and we can have the same equation for ip by complex conjugating of equation 
(9). We take the ansatz: 



ij = V(r), At = 0(r), Aa = 0,a = r, 9, 



(11) 
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It is then suitable to take the phase of ip to be constant. All other fields are set to be zero. 
Under this ansatz, the equations of motion simplify to: 

^7-2(^2-7 j^')' + + N-^(f>'^i: = (12) 

7-2^^2+7^/y _ 2(j)^\^y-^ = (13) 



r 



where a prime denotes the derivative with respect to r, and we have to notify that if 7 = 0, 
these equations will reduce to the ones in [l|, 0, 29 1. We define a mass parameter as: 



m^L^ = -2 

The field equations (12), (13) can be written as the next set: 

(f)" + (2 + 7)r2T'- V - 20VV2T/-1 = (15) 

If /3 = 1, 7 = we recover again the results of |l|, 0, We must note an important fact 
about the limiting process to achieve the Lu et al solution given in js^. The limiting process 
7 — )■ is valid for both different values of the A = 1, 3 > 1. The Lu et al solution recovers 
both of these values, although we observe from the form of the lapse function that these 
values lead to the same metric functions . 

Examining these fields equations at the horizon and assuming that the scalar field must 
be regular on the horizon, we can observe that we have the next set of the auxiliary boundary 
conditions: 

= ^^r, (16) 
Jh 



(17) 
in which th is the horizon radius of the black hole, i.e. the largest root of /(r) = 0. 

A. Solving the general equations in the asymptotic region 

In the vicinity of the black hole, Eqs (14), (15) can be solved by making a change of 



variable, r vh, and setting the radius of AdS"^ to be L = 1 [30]. In [30| also the case 
= was discussed both via numerical and semi analytical methods. In this manuscript 
we limited ourselves only to this special case = 0. We can easily guess their behavior in 



8 



the large r limit. In order to find the asymptotic behavior of the field we must determine 
when in the IR region A > 1, the exponent 6 is positive or negative. There are two different 
kinds of the exponent 6 which we denote them by 5+,5„. We mention here that for a 
sufficient large value of the A the value of the exponent S- remains below 2. Thus for all 
values of the A > 1, we have the next limiting values: 



lim (5+) 




+ 00 


(18) 


lim 




1 

2 


(19) 




< 


2 


(20) 


2 < (5+ 


< 


oo 


(21) 


1 < 7+ 


< 


oo 


(22) 


1 

-2 < 


< 


1 


(23) 



B. Approximation techniques 



According to the method discussed in 3l| we must find the approximate solutions near the 
horizon, then generalize it to the asymptotic AdS region and smoothly match the solutions 
at an intermediate point. By introducing a new radial-like coordinate as: 

e = - (24) 
r 

we can rewrite the equations (14), (15) in terms of the new coordinate ^^: 

^ f . 27+4^.-27-4 



+ (7 + 7)^ + ( at, ^')^ = (25) 



+ (| - 4^(2 + 7)r^-^^)0 - 2^H^+\-2^-V- V = (26) 

where a dot now denotes 4: and we observe that for the interval out of the horizon this 
coordinate smoothly covers all points of the strip: 

th <r <oo, 0<^<1 (27) 

The boundary conditions (16) and (17) in the massless limit with the regularity at the 
horizon ^ = 1 become: 

0(1) = 0, V'(l) = (28) 



We limited ourselves to a massless case = 
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With this change of the variable the equations (14) and (15) convert to the next set (16) 
and (17), which must be solve near horizon i.e ^ = 1 with auxiliary boundary conditions 
(28). Our main goal is to find the coefficients and powers in (25), (26) and also matching 
these two solution in an intermediate point. 



C. Solutions near the horizon:^ = 1 

We can expand ip{r) and 0(r) in a Taylor series near the horizon as: 

1 



0(0 = 0(1) - 0(1)(1 - + 2'^(1)(1 - 0' + - (29) 

1 v. 



^'(o = V'(i) - - + ^V'(i)(i - 0' + - (30) 

According to the equation (28), for a massless scalar field, we have -0(1) = and (f)(1) — 0, 
and without loss of generality we take 0(1) < 0, ■0(1) > to have 0(1) and ^0(1) positive. 
Expanding (26) near ^ = 1 gives: 

0(1) = (^rr^ + rl\2 + 7) - 2)0(1) (31) 

Thus, we get the approximate solution: 

m = 0(1)(-(1 -0 + 2(1 - 0\ ^^-^^^ + 4"(2 + 7) - 2)) (32) 
Similarly, from (25), the 2'nd order coefficients of can be calculated as: 

^(1) - -l^^dxlj)^ (33) 

where we used Hopital rule at the second term, therefore an approximate solution near the 
horizon is: 

^(0 = ^(i)(i-^(|^f)(i-e)^ (34) 

D. Solutions in the asymptotic AdS region 

In the asymptotic AdS region ^ = 0, the solutions are: 

ii^ ^ D^e^ + D^e- (35) 

(j)^ li-q^ (36) 
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where fi is the chemical potential and q is the charge density on the boundary^ . At the 
boundary of a (2+l)-dimensional field theory, fi is of mass dimension one and q = p/th 



is of mass dimension two. From the boundary be 



fly 



laviors, we can read off the expectation 
3211 ■ we know that , both of these falloffs 



value of operator O dual to the field. From ^ ^ 

are normalizable, and in order to keep the theory stable [l|, we should impose the following 
equations: 

= 0, < d_ >= V2D_ (37) 
D_ = 0, < 6+ >= V2D+ (38) 

where the factor \/2 is a convenient normalization Ij . The index i in Di represents the 
scaling dimension Xq of its dual operator < Oj >, i.e. Aq. = i. Note that these are not 
entirely free parameters, as there is a scaling degree of freedom in the equations of motion. 
As in [1], we impose that p is fixed, which determines the scale of this system. For ip, both 
of these falloffs are normalizable, so we can impose the condition either D_ or vanish. 
We take D_ = 0, for simplicity. Now we must find the solutions of the equations (25) and 
(26) with the boundary conditions mentioned above. Since the dimension of temperature 
T is of mass dimension one, the ratio T^/p is dimensionless. Therefore increasing p, while 
T is fixed, is equivalent to decrease T while p is fixed. We must show that when p > p^, 
the operator condensate will appear; this means when T < Tc , there will be an operator 
condensation, that is to say the superconducting phase occurs. 

We limited ourselves only to the case 5+ > 2, 7+ > 1. Remembering for a general 2'nd order 
differential equation, we can write (25) in the following self-adjoint form: 

+ P(a;)^ + Q(a;)* = (39) 

The change of the variable \E'(x) = e~^/^-/^'^*^^)'^^H(x) converts it to the next Schrodinger like 
equation: 

H(x) + (-1/2P - 1/4P2 + Q)E{x) = (40) 

For (25) we have: 



Our compendium follows what mentioned in the Gregory et.al work [3 
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In AdS asymptotic region with the metric function f ^ —a^ ^, the field equation (40) is 
converted to the: 

em + vm) = o, 77 = i/4(i -{^-5- ir) = ~ (42) 

This is a standard Euler-Cauchy equation which has the following exact solution: 

s(e) = s+r+ + s_r-, ^± = |-^ (43) 

m = D+e + D- (44) 

m -y^-pi (45) 

The new set of coefficients D± are some functions of the S±, a, .... 



IV. MATCHING AND PHASE TRANSITION 

Now we will match the solutions (32), (34), and (44), (45) at Allowing to be 
arbitrary does not change qualitative features of the analytic approximation, and more 
importantly, it does not give a big difference in numerical values; therefore for simplicity in 
demonstrating our argument we will take = 1/2. In order to connect our two asymptotic 
solutions smoothly, we require continuity in our fields and their first derivatives at the 
crossing point = 1/2, therefore following four conditions should be satisfied^: 

^_^ = _6(_5+l(_^^+47(2 + ^))) (48) 

/' = ('(2-i(-^^^|^+r|'(2 + 7))) (49) 

after setting D_ = 0, we obtain from equations (46) and (47): 

= 2« = (50) 



(51) 



b ^ bT (52) 



We have set ^(1) = a and -0(1) = b, {a, 6 > 0) for clarity, /(I) = -^/i^+i 
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FIG. 1: A plot of the critical temperature as a function of p and 7 varying in the range of 
-0.5 < 7 < 1 



where {h := ^^^ ) and also from equations (48) and (49) we have: 

2 167rT ^ p 3,,^ /i27(2 + 7) 

^=m^^^^~2-t^' 



bl3h'r+3 

where ( h := Vh) and then we conclude that: 



4 



b = ip, — 3p 



and we can define the critical point, Tc as: 



26(1 



4 



(53) 
(54) 

(55) 

(56) 



Figure (1) shows the the dependence of as a function of p and 7. As we see when 
p = for different values of 7 the value of Tc is equal to zero, and in the case 7 = there is 
a linear dependency of T^. with respect to the varying parameter p. This is also mentioned 
in the figure (2) . As we see in the figure (2) when p = the magnitude of is equal to zero 
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FIG. 2: A plot of the Tc as a function of p ( 7 = ). In this plot 7 = as we change p between 
and 10. As we see there is linear dependency with respect to parameter p. 

and when p goes higher the Tc also goes higher with Hnear dependency. In the figure (3) we 
show the dependency of Tc with respect to 7 in the range of —0.5 < 7 < 1, when p is fixed 
(for example in that case p = 10). With increasing of 7, the values of Tc also increase but 
not linearity. 

Noting that in order to remain the temperature Tc positive, we must have h'^ < ^2+^ 
, and according to the equation(7) we can conclude that(| < h < , and it could be 
reasonable to choose h = 1 . Near the critical temperature the AdS/CFT dictionary gives 
the relation below : 



< 0+ >= V2D^ = 2V2a = 4./ - ^Y^' (57) 

We observe that < 0+ > is zero at T = Tc, the critical point, and condensation occurs 
for T < Tc- The continuity of the transition can be checked by computing the free energy[l|. 
We also see a behavior < 0+ >oc [Tc — T)^/^ which is a typical mean field theory result for 
a second order phase transition j33|. 

Figure (4) shows (0+) as a function of temperature normalizing by Tc for a variety of 
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FIG. 3: A plot of the Tc as a function of 7 ( /) = 10 ). This plot shows how by varying 7 the 
values of Tc change when p is fixed ( for this case p = 10 ). 

values of p and (3. Each line in the plot forms the characteristic curve of (0+) condensing 
at some critical temperature. For simplicity we chose five values of p and /3 to display the 
features of the system and showing how varying /3 and p effect the height of (0+). In this 
figure according to the equation (7), in order to have positive Unruh temperature we must 
require A^y < 0. 

As we see in figure (4) increasing /3 reduces the value of We also see that the 

condensation appears when T = Tc. 

Figure (5) shows that the effect of increasing p is to increase the height of these graphs 
( (0+) ), in similar way mentioned in figure (4), the condensation happens at T = T^. 

V. CONDUCTIVITY 

In order to compute the electric conductivity in dual CFT , we must solve the Maxwell 
equation for the fluctuations of the vector potential Ax{r, t), located in the bulk. We assume 
that the time dependence of the field is e"*"^* and then the field equation of this component 
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p=i,p.i 

p.1,p.1.5 
p=1.p.2 
p.1,p=2.5 

p.l,i;.3 




T/To 



FIG. 4: The plot of the condensation as a function of y' ^ selection of values of p and (3. In 
this plot the value of p is fixed ( in this case is equal to one) and the value of /3 from top to down 
is equal to 1, 1.5, 2, 2.5, 3. As we see the height of (0_|_) is decreasing as the values of /3 increase, 
and the condensation occurs when T < Tr. 



reads as: 



which is what mentioned in the papers 



33 



T 



A = 



(58) 



35| in the special case k = and a 



where from the metric ansatz we have concluded that 



I3r- 



The causal behavior is 



obtained with imposing an ingoing wave boundary condition at the horizon 
asymptotic behavior of the Maxwell field at large distance is 



. The desired 



A 



(1) 



(59) 



According to the AdS/CFT dictionary, the dual source and expectation value for the current 
are given by 



(60) 
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FIG. 5: The plot of the condensation as a function of y' ^ selection of values of p and /?. In 
this plot the value of f3 is fixed ( in this case is equal to one) and the value of p from down to top 
is equal to 1, 1.5, 2, 2.5, 3. As we see the height of (0+) is increasing as the values of /3 increase, 
and the condensation occurs when T < T^. 

Now using Ohm's law we can obtain the conductivity as 

.M^^ (61) 

Thus we must solve (58) numerically and obtain the imaginary part of the conductivity cr(c(;) 
for a set of parameters Aw = —2, P = —1,6 = —1. There is a delta function at w = which 
appears as T < Tc, and from the Kramers-Kronig relation we can see that the real part 
of the conductivity contains a delta function and the imaginary part has a simple pole at 
a; = 0. Thus the superfluid density is of the delta function [l| 

Re{a{uj)) ~ 7rS{uj) (62) 

The figure (6) shows the behavior of the real part of conductivity as a function of frequency 
per temperature for different values of the HL parameter a. 

We can solve (58) for operator O2 for the former set of the parameters. The result graph 



17 



1.2 
1.0 
0.8 

Re(o-) 0.6 
0.4 
0.2 
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I 

FIG. 6: The real part of conductivity Re{(T{uj)) as a function of the frequency ^ for Oi operator 

1.2 
1.0 
0.8 

Re(o-) 0.6 
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0.2 
0.0 

50 100 150 200 

&) 

T 

FIG. 7: The real part of conductivity Re{(T{uj)) as a function of the frequency ^ for O2 operator 
has been shown in the figure (7). 

VI. CONCLUSION 

In the present work, we have built a holographic model for a non-relativistic system 
showing superconductivity. We have used a black hole background which comes from the 
Horava-Lifshitz gravity, and we have studied analytically, holographic superconductors in 
this new kind of the asymptotic AdS solutions. We also have analytically solved the system 




50 100 150 200 



T 
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in the probe limits, near horizon and asymptotic region. We have found that there is also 
a critical temperature like the relativistic case, below which a charged condensation field 
appears by a second order phase transition, and also we have found out below a critical 
temperature Tc, the condensation field appears and obtains finite value. We can conclude 
that as the condensation field becomes heavier, the transition happens more observable. 
Also the conductivity has been computed and the variation of the critical temperature and 
conductivity with respect to the parameters of the metric function have been shown. We 
numerically obtain the conductivity as a function of the frequency for a wide range of the 
parameters. We show that the Gauss-Bonnet theory in five dimension and Horava-Lifshitz 
theory in critical exponent z = 3 and in four dimension share some similar features. 
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